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Abstract Currently, both the g factor measurement of the muon as well as the Lamb shift 2S-2P 
measurement in muonic hydrogen are in disagreement with theory. Here, we investigate possible 
theoretical explanations, including proton structure effects and small modifications of the vacuum 
polarization potential. In particular, we investigate a conceivable small modification of the spec- 
tral function of vacuum polarization in between the electron and muon energy scales due to a 
virtual millicharged particle and due to an unstable vector boson originating from a hidden sector 
of an extended standard model. We find that a virtual millicharged particle which could explain 
the muonic Lamb shift discrepancy alters theoretical predictions for the muon anomalous mag- 
netic moment by many standard deviations and therefore is in conflict with experiment. Also, we 
find no parameterizations of an unstable virtual vector boson which could simultaneously explain 
both "muonic" discrepancies without significantly altering theoretical predictions for electronic 
hydrogen, where theory and experiment currently are in excellent agreement. A process-dependent 
correction involving electron screening is evaluated to have the right sign and order-of-magnitudc 
to explain the observed effect in muonic hydrogen. Additional experimental evidence from light 
muonic atoms and ions is needed in order to reach further clarification. 
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Section 1: Introduction 



1 Introduction 

Recently, two experiments involving quantum electrodynamics (QED) effects have been in disagreement 
with theory. The muon anomalous magnetic moment a M = (g^ — 2)/2 has been measured [IH3] as 

a M cxp = 11659208.0(6.3) x 10~ 10 (1.1) 

in 3.4<t disagreement with some of the latest theoretical analyses [4] 

a p th = 11659180.4(5.1) x 10~ 10 . (1.2) 

The original aim of the recent muonic hydrogen Lamb shift experiment [5] was the determination of the 
proton radius. When QED theory is assumed to be correct, then the value 



2) p = 0.84184(67) fm (1.3) 

is inferred for the root-mean-square proton charge radius from a comparison of theory and experiment 
for the transition 2S\/2{F — 1) 2P 3 / 2 (F = 2) in muonic hydrogen (/xH). This value of the proton 
radius is in disagreement with the value obtained in the same way mainly from hydrogen and deuterium 
spectroscopy [6] , which is the basis of the CODATA value [7] , 

r p = 0.8768(69) fm. (1.4) 

The most recent and accurate measurement of the proton radius from electron scattering jS], yields a 
value of 

r p = 0.879(8) fm, (1.5) 

when the statistical and systematic uncertainties given in Ref. [8] are added quadratically. The two 
values (|1.4[) and (|1.5|) are in excellent mutual agreement but differ from the muonic hydrogen value (|1.3|) 
by 5.0 standard deviations. Consequently, it may be permissible to invert the argument, and to evaluate 
current QED theory for the muonic hydrogen transition (as summarized in the supplementary material 
published with Ref. [5]) with the CODATA value (|1.4p . Using the theoretical expression given in Ref. [5], 

_ ( r 2 r 3 \ 

Eth = 209.9779(49) - 5.2262-^ + 0.0347-^ mcV, (1.6) 
V fm fm / 

one obtains, using the CODATA proton radius given in Eq. (jl.4p . a theoretical prediction of 

E th = 205.984(63) fin, (1.7) 

which is 5.0cr away from the experimental value of 

E cxp = 206.2949(32) fm , (1.8) 

reported in Ref. [S]. The recent theory update [S] shifts theoretical predictions only minimally, to 

E th = ^209.9974(48) - 5.2262^-^ meV . (1.9) 

Using the CODATA proton radius given in Eq. (|1.4p . one then obtains the theoretical prediction of 

E th = 205.980(63) fm, (1.10) 

in excellent agreement with (11.71) . but in significant disagreement with the experimental result given in 
Eq. (HU). 
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Section 2: Historical Perspective on Discrepancies in Muonic Systems 



Both of the most recent QED experiments involving muons [3][5] are in disagreement with theory. The 
discrepancies have the "same sign" and read 

Sa„ = a M oxp - a M th = 2.76(81) x 1CT 10 , (1.11a) 
SE = E CKp - E th = 0.316(63) meV . (1.11b) 

We here proceed as follows. First, in Sec. [2 we present a historical perspective on discrepancies in muonic 
bound systems observed in the past, and on their eventual resolution. In view of the current discrepancy, 
such a historical perspective may be useful. 

Possible theoretical explanations for the current discrepancy can mainly be divided into two categories: 
proton structure effects and modifications of the vacuum polarization charge density. Hypothetical proton 
structure effects are discussed in Sec. [3j and modifications of the vacuum polarization charge density in 
Sec. |U Section [3] is divided into two subsections, the first of which deals with a conceivable "dip" in 
the proton form factor slope in the momentum transfer range studied in muonic hydrogen spectroscopy, 
and the second deals with a conceivable, anomalously large contribution from the inelastic part of the 
two-photon exchange diagram. Section U is divided into three parts. These deal with a conceivable non- 
pcrturbative correction to the vacuum polarization potential, with the contribution of a millichargcd 
particle that modifies the vacuum polarization loop, and with the contribution of a conceivable virtual, 
unstable vector boson that modifies the vacuum polarization potential. We anticipate here that none 
of these considerations will lead to a definitive candidate for an explanation of the discrepancy. Still, a 
compilation of a number of possible explanations appears to be useful in the current situation. Some of 
the discussed explanations are relevant only to muonic hydrogen, which is the main subject of the current 
article, others may be relevant for both observed discrepancies. Possible electron screening corrections 
are discussed in Sec. 03 Conclusions are reserved for Sec. [6] Natural units with H = c = 6q = 1 are used 
throughout the paper. 



2 Historical Perspective on Discrepancies in Muonic Systems 

Let us briefly comment on the size of the disagreement in muonic hydrogen [see Eqs. (|1.10[) and (jl.8l) ]. 
Current predictions are based on the calculations reported in Rcfs. [TUHT7] and represent the result of 
independent groups. Important contributions originally calculated in Rcfs. llOUllj have been verified in 
Refs. fl3l[14] . Higher-order vacuum polarization effects have been given special attention in Refs. [T3l[T4]. 
The theory used in the evaluation of the experiment [5] has been compiled at Laboratoire Kastlcr-Brosscl 
in Paris. The disagreement of theory and experiment is on the level of ^ of the two-loop vacuum polar- 
ization correction which amounts to 1.508 meV and which was calculated first by Kallen and Sabry [18] . 
then recalculated in Ref. [19]. A clear exposition is given in volume III of Ref. [20] (the result has later 
been generalized to non-Abelian gauge theories, see Refs. [21T - I23] ). 

One may point out that the discrepancy (|l.llb[) amounts to (roughly) 1.5 parts per thousand of the 
total vacuum polarization effect in muonic hydrogen. By contrast, in a previous experiment |24] involving 
muonic 3c? 2p transitions in 24 Mg and 28 Si, the vacuum polarization effect has already been verified 
to 1.0 parts per thousand (a relative accuracy of 950 x 10~ 6 is quoted in Ref. [24]) and thus, to better 
precision than the current disagreement. If there were any fundamental reason for a deviation of theory 
and experiment on this level, then one might wonder if the effect (whichever it is) might have been visible 
in the experiment reported in Ref. |24j . However, the muonic transition reported in Ref. |24j suffers from 
an uncertainty due to electron screening, and also, as it involves non-S* states, the overlap of the muonic 
wave functions with the nucleus are not as pronounced as for S states. So, the quoted experiment |24j 
is not sensitive to higher-order nuclear structure effects, and it also probes the vacuum polarization at a 
different energy scale as compared to S states which are much closer to the nucleus. 

In the 1970s, experiments involving muonic transitions were found to be in disagreement with theory |25j . 
Part of the discrepancies were addressed after a sign error in the calculation of the two-loop vacuum 
polarization correction [26] was eliminated [27]. An elucidating discussion of the status reached in 1978 
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Section 3: Proton Structure Effects 



is given in Ref. |28j . Further clarification was reached when a standard 7-ray spectrometer used in the 
experiments was recalibrated |29j . A few remarkable experiments later found full agreement of theory 
and experiment in muonic systems (e.g., Refs. ) - In Ref. [31], the nuclear radii of some carbon, 

nitrogen and oxygen isotopes are determined by analyzing muonic transitions, and the resulting radii are 
in agreement with electron scattering radii for the investigated nuclei to better than 5 %. Later, the radius 
of 12 C has been updated in Refs. |32l33j . converging at a value of = 2.478(9) fm, in excellent agreement 
with the value from muonic x-ray studies. After this finding, a conceivable non-universal coupling of the 
electron versus muon to the nucleus was discarded, and muonic spectroscopy meanwhile is an established 
tool for the determination of nuclear radii |34j . 

Indeed, the current disagreement of the proton radius derived from muonic versus electronic hydrogen 
radius amounts to about 4% and is a large discrepancy (roughly, electron scattering and hydrogen 
spectroscopy gives a result of r p w 0.88 fm, whereas the recent muonic Lamb shift experiment yields 
r p m 0.84 fm). Therefore, it is permissible to speculate about nonperturbative effects and new physics 
effects for both discrepancies (muon g factor and muonic hydrogen Lamb shift). Among the two effects, 
the hydrogen experiment is perhaps the most interesting, (i) because its theory, on the level of the 
discrepancy, is given by only few, simple bound-state QED effects not exceeding the two-loop level (see 
Ref. [9]), and (ii) because new experiments in related systems are planned which may or may not confirm 
the observed discrepancy [35]. Here, we thus investigate a few of these possible explanations. 



3 Proton Structure Effects 



3.1 Form Factor 



The proton mean-square charge radius is defined in terms of the slope of the Sachs Ge form factor of the 
proton, 

dG E (q 2 ) a dG E {Q 2 ) 



(r 2 ) = 6 



6 dQ 2 



(3.1) 

Q 2 =0 



q 2 =0 

where Q 2 = —q 2 is the space-like momentum transfer. Different ranges of the momentum transfer are 
relevant for the calculation of the slope in different experiments. The proton radius from electronic 
hydrogen is determined from exchanged Coulomb photons with momentum transfers in the region 



c) 2 = ( 3.7 x lO- 6 ^^ . (3.2) 



For muonic hydrogen, the atomic momentum is in the range of 

Q 2 ~ {am^cf = ^7.7 x 10~ 4 ■ (3-3) 

One may point out that this is just below the electron-positron pair production threshold, 

, GcV\ 2 
\^ — 1 n v 1 n— J 1 



(2m e cY = 1.0 x 10~ d . (3.4) 



c 



The momentum transfer range probed in the recent electron scattering experiment [8] is larger, but not 
excessively larger, 



, 2 GeVV ^2 ( GeV 



c 



6.3x10^—1 <Cr<(l— ) . (3.5) 



2 



The slopes of the proton form factor determined from the electron scattering data and from electronic 
hydrogen spectroscopy are in excellent mutual agreement. The momentum transfer range for muonic 
hydrogen spectroscopy lies in between these two ranges. Consequently, it would be somewhat surprising 
if the proton form factor slope had a "dip" in this range that would explain the discrepancy for muonic 
hydrogen. Still, without a direct scattering measurement in this momentum transfer range, this possibility 
cannot be fully excluded at present. 
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Section 4: Vacuum Polarization 



3.2 Proton Polarizability 

Generically, the proton polarizability can be related to the resonances of the proton (its excitation spec- 
trum) via dispersion relations [see Eqs. (29) and (30) of Ref. pT]], and related to the inclusive reaction 
e+p — > e' + X. This is an accepted procedure for all nuclei, also for heavier nuclei [3T]. For the proton, one 
would intuitively assume that the bulk of the contribution is from the A(1232) resonance, which has been 
measured well. Yet, the data obtained in the literature for the proton polarizability contribution scatter, 
and in Ref. [5], the contribution is currently estimated as +0.015(4) meV based on the scatter of values 
obtained from different theorists [IT][36l|37]. The proton polarizability correction to the 2P-2S Lamb 
shift has been calculated as 0.0012 me V in Ref. [TT], which is an order to magnitude smaller than the 
discrepancy SE. Other authors |31)Il3"7] confirm the magnitude of the result and give values of 0.0015 meV 
for the proton polarizability correction to the Lamb shift in muonic hydrogen. It would thus be helpful 
to reevaluate the effect, and to obtain more accurate estimates, even if the current uncertainty estimate 
of ±0.004 meV is numerically tiny as compared to the discrepancy of SE w 0.31 meV. 

In addition, one may point out that in the past, nuclear radii inferred from muonic transitions and from 
electron scattering have agreed to better than 5% (see the discussion in Sec. [5]). So, if an anomalously 
large proton polarizability contribution were found by a reanalysis, then one might have to revisit this 
effect also for other bound systems and in the more general context of the validity of the nuclear charge 
radius determination from muonic transitions |34U38) . 

4 Vacuum Polarization 

4.1 Nonperturbative Vacuum Polarization 

Let us briefly review why a nonperturbative vacuum polarization effect might have been considered as an 
explanation of the discrepancy in muonic hydrogen. The spectrum of muonic hydrogen is influenced by 
electronic vacuum polarization effects. The muon is heavier than the electron by a factor of m M /m e w 207, 
and the reduced mass of muonic hydrogen is roughly equal to the muon mass. The effective Bohr radius 
in muonic hydrogen is l/{ams) = 284.748 fm, which is smaller than the reduced Compton wavelength of 
the electron, l/m e = 386.159 fm. The bound muon thus enters the electronic vacuum polarization charge 
cloud of the proton. The electronic vacuum polarization shift in muonic hydrogen is of order a 3 m/j, 
where is the reduced mass, and thus more pronounced than the electronic vacuum polarization shift 
in electronic hydrogen, where the vacuum polarization contribution to the Lamb shift is of order a 5 mR. 
The large vacuum polarization shift is also responsible for the fact that the 2S level in muonic hydrogen 
is energetically lower than the 2Pi/ 2 state (in contrast to electronic hydrogen, where the situation is 
opposite). 

Superficially, the vacuum polarization effects converge very well in terms of the QED loop expansion. The 
one-loop effect gives a contribution of 205.0074 meV to the 2P-2S Lamb shift, in first-order perturbation 
theory, while the second-order effect adds 0.1509 meV. The two-loop (Kallen-Sabry) shift is 1.5081 meV, 
followed by the higher-order Wichmann-Kroll term of —0.00103 meV. Eventually, of course, the expansion 
will diverge according to a famous argument put forward by Dyson |39j , but the expected nonperturbative 
effect would be of order exp(— 1/a) and thus completely negligible for the muonic hydrogen experiment. 

However, the superficial convergence still does not exclude the presence of a much larger nonperturbative 
correction to the vacuum polarization if the local convergence of the loop expansion of the higher-order 
vacuum polarization potentials breaks down in close vicinity of the proton, i.e., if the higher-order (Kallen- 
Sabry and Wichmann-Kroll) terms are more singular than the one-loop Uehling term for r — > 0. In 
that case, a nonperturbative correction to the vacuum polarization potential might have led to a highly 
nonlinear, nonperturbative correction and one might have had to solve the Schrodinger equation using 
the full nonperturbative potential near the origin. In that case, lattice methods would probably have had 
to be invoked in order to calculate the full vacuum polarization potential for r — > 0. 

The question whether this more elaborate calculation is necessary, can only be answered by a concrete 
calculation of the leading asymptotics of the Uehling, Kallen-Sabry, and Wichmann-Kroll potentials 
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Section 4: Vacuum Polarization 




Figure 1: Vertex correction to the muon anomalous magnetic moment 
with an additional vacuum-polarization insertion in the photon line (left: 
electron-positron loop, right: non-QED virtual process). 



for r — > 0. One finds, in agreement with Rcf. [57], for the leading asymptotics of the one-loop Uchling 
potential V vp (r), 

Mr)~!^]n(m e r), r^O, (4.1) 

for the Kallen-Sabry potential, 

4 a 3 , 

Vks r) ~-— In 2 (m e r), r^O, (4.2) 
9 irr 

and for the Wichmann-Kroll potential, 

7 WK (r)~^ (-§C(3) + |tt 2 - , r^O. (4.3) 

These potentials have to be compared to the Coulomb potential 

V(r) = --. (4.4) 
r 

By inspection of these formulas, we conclude that the higher-order vacuum-polarization potentials are of 
the same order-of-magnitude as the Coulomb potential for distances shorter than 

r ^exp(-l/q) =12xl0 _ 57£m! 

which is the length scale of the Landau pole. This length scale is not sufficient to induce to any conceivably 
large nonpcrturbativc effects. 



4.2 Virtual Millicharged Particles 

One of the most interesting possibilities for an explanation of both "muonic QED discrepancies" observed 
at present would be due to the contribution of a virtual millicharged particle. A millicharged particle was 
invoked as a possible explanation for the observed (later retracted) optical rotation [40H42] of linearly 
polarized laser light by a magnetic field. If the photon initiates pair production of light charged fermions 
with masses below the electron mass and charge on the order of q = ee with e <C 1, then the initial 
observation made in Ref. [ID] could be explained (see Refs. [4*3"H4"5] ). The non-integer charge does not 
contradict charge quantization if the millicharged particles are generated from a "hidden" sector of the 
standard model via the Stueckelberg mechanism [46][47] . Such millicharged particles have been searched 
in devoted experiments at SLAC [35] . Some of these experiments are sensitive only to stable millicharged 
particles, because they depend on obtaining a signal from particle detectors, as pointed out in Ref. |49j . 
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Section 4: Vacuum Polarization 




Figure 2: Vacuum polarization correction to the muonic hydrogen Lamb 
shift (left: electron-positron loop, right: the blob denotes a correction due 
to a non-QED virtual process, e.g., hadronic vacuum polarization, p(770) 
meson pole, or due to a hypothetical low-energy vector meson). 



Several bounds have been derived regarding the mass and charge of such particles, which have otherwise 
been quoted as a candidate for dark matter (see Refs. [5DH55] ). 

The muon anomalous magnetic moment discrepancy and the Lamb shift discrepancy have the same sign, 
i.e., the experimental result is larger than the theoretical prediction. An additional virtual excitation of a 
quantum field (a virtual particle) would naturally be assumed to enhance both effects. The muon anoma- 
lous magnetic moment is numerically small, the correction induced by a hypothetical virtual particle is 
a two- loop effect (see Fig. [T]), whereas for the muonic hydrogen Lamb shift, the conceivable contribution 
of a millicharged particle only is a one-loop correction (see Fig. [2]). Therefore, it is indicated to map out 
possible parameter ranges for the hypothetical millicharged particle. In muonic experiments, one is very 
sensitive to the mass range m e -C mj/ <C for hypothetical virtual particles. If the virtual particle is 
in this mass range, then the effect on the muon anomalous magnetic moment, and on the muonic Lamb 
shift is enhanced because tum <C tu^, but suppressed for electronic systems such as ordinary hydrogen 
because m e <C Tom- These models are simple-minded, straightforward ansatzes that "suggest themselves" 
because of the pertinent mass region. We thus restrict the discussion to conceivable millicharged particles 
and do not consider supersymmetric graphs in which, e.g, the muon might turn into a virtual smuino, 
emitting a charged higgsino or wino. We also do not consider hypothetical corrections from axion elec- 
trodynamics [55] , 

Furthermore, since m e <C tum by assumption, the modification of the vacuum polarization due to the 
millicharged particle, for electronic hydrogen, can be absorbed into a Dirac 5 potential acting on the 
electronic hydrogen wave functions. Its functional form is therefore indistinguishable from the nuclear 
finite size effect for electronic hydrogen and could be "absorbed into" a modification of the proton radius 
inferred from electronic hydrogen spectroscopy without any further observable consequences for atomic 
transitions in electronic hydrogen. For muonic hydrogen, however, since mj/ <C m^, the hypothetical 
millicharged particle leads to an enhanced energy correction, different from a Dirac S. The mass range 
m e <C mm *C to m therefore is the primary parameter range probed by muonic QED experiments. For 
to a/ 3> to m ^> TO e , the effect of the millicharged particle amounts to a Dirac 5 function and is thus 
indistinguishable from the contribution of the nuclear size effect for both muonic as well as electronic 
hydrogen. 

We thus proceed as follows. First, we find a convenient parameterization of the expected modification of 
the spectral function of vacuum polarization, as a function of the charge and mass of the millicharged 
particle. Then, evaluate the shift of the anomalous magnetic moment of the muon and of the Lamb 
shift due to the millicharged particle (these are both proportional to the square of the charge) , and relate 
them to the observed discrepancies, as a function of the mass of the assumed millicharged particle. Again, 
forming the ratio of these relative shifts, we investigate if there is a parameter range for which both the 
anomalous magnetic moment of the muon and the discrepancy observed in muonic hydrogen could be 
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Section 4: Vacuum Polarization 



G(m M ) 




20 40 60 80 m t 



Figure 3: (Color online.) In the range 10 m e < tum < 100 m e , the function 
G(uim) is a lot larger than unity, as shown in the plot. If a hypothetical 
millichargcd particle in the given mass range were responsible for the 
discrepancy observed in muonic hydrogen, then the same particle would 
lead to complete disagreement for the anomalous magnetic moment of the 
muon. 



explained by the virtual particle. 

First, let us find a convenient parameterization for the spectral density of vacuum polarization. As is well 
known, the effect of electronic vacuum polarization on the photon propagator can be described by the 
replacement 

1 a f°° dt , . 1 

2 , • T~ / TPeW-2 — TT^' 
<T + ie Sir J 4m 2 t q A - t + ie 

/ Ami ( 2m 2 \ 

^(*) = Y 1 "~r ( 1+ ~rJ ' (4 - 6) 

in the photon propagator, where we refer to p e (t) as the spectral density. The corresponding vacuum 
polarization potential (Uchling potential) induced by electronic vacuum polarization is 

K P (r) = -— / Pe (t). (4.7) 

J7T J im 2 t r 

The spectral density is zero at threshold t = and quickly approaches the asymptotic value p e (t) — > 1 
for t — > oo. Although the millicharged particle has been assumed to be of spin 1/2, we emphasize here 
that similar threshold behavior can be expected regardless of the spin of the millicharged particle [571158] . 

A millicharged particle modifies the spectral density of vacuum polarization according to p e (t) — > p e {t) + 
5p(t). For a spin-1/2 millicharged particle, of charge q = ee and mass m,M, we can approximate this 
modification as 

5p(t) w e 2 6(f - Am 2 M ) (4.8) 

where is the step function. We intend to compare changes in the anomalous magnetic moment of the 
muon and in the Lamb shift induced by the millicharged particle. In line with intuition, we here need a 
positive spectral function Sp(t), increasing the Lamb shift for muonic systems and increasing the muon 
g factor. This is the right sign because the theoretical prediction for the muon g factor as well as the 



8 



Section 4: Vacuum Polarization 




Figure 4: (Color online.) Both the muon anomalous magnetic moment 
discrepancy as well as the muonic hydrogen Lamb shift discrepancy can 
be explained by a millichargcd particle of mass hum = 0.221 m e and charge 
q = ±0.0179e, as shown in the graph. Indeed, one finds G(0.221 m e ) = 1. 
However, in the indicated mass range, the correction to the electronic 
hydrogen Lamb shift induced by the millicharged particle becomes so 
large that it leads to an inconsistent, sizeable shift of the the proton radius 
inferred from the hydrogen Lamb shift. See text for further explanations. 



theoretical prediction for the muonic helium Lamb shift are lower than the corresponding experimental 
results. 



The correction to the anomalous magnetic moment due to electronic vacuum polarization is |59U60| 

a 2 f°° dt , N , , . , 

= 7 2 / T fa(t) P e(t), (4.9) 
b7T Ural * 

where for the muon 

We have checked that if one replaces in this expression — > m e and integrates over t, then one obtains 
the known contribution [61] . 

Sa e = — - T , (4.11) 

to the electron anomalous magnetic moment a e , due to the diagram on the left in Fig. [T] 

Although the integral representation of f a {t) is compact, the analytic result requires us to differentiate 
two cases, depending on whether t < 4m^ or t > 4m^. For < t < 4m^, with r = tj (4m^), one finds 

f a {t) = 1 -St- 8r(l - 2r) ln(4r) -4(1 -8r + 8r 2 ) -^j arctan (\f^T^ J - (4.12a) 

whereas above the muon threshold, for t > 4m 2 , , with x = [1— \ 1 — Am 2 It) / ( 1 + . 1 — Am 2 It], the 



result for f a {t) is 



(l - y/l - I (l + y/l-Aml/t), 



9 / , , / ln(l + x) — x + i x 2 l + x 
f a (t) = x 2 2 - x 2 + 2 1 + x 2 1 + a; 2 -i ^ 2 - + 2 x 2 In a; • 

x z 1 — x 
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Section 4: Vacuum Polarization 



Consequently, the contribution to the muon anomalous magnetic moment due to the hypothetical mil- 
licharged particle, divided by the observed discrepancy 5a given in Eq. (|l.lla[) . is 

Xa = jT°°Y »?,(«) (4.13) 

where 

Here, Sp{t) is the vacuum polarization spectral density given in Eq. (|4.8|) due to the millicharged particle. 

For the muonic hydrogen Lamb shift, the situation is as follows. The one-loop electronic vacuum po- 
larization shifts the 2S level downward by — 219.6 mcV, whereas the 2P is shifted downward by only 
— 14.6meV. This is because of the enhanced probability density of the 25 state as compared to 2P, 
near the nucleus. The total effect on the Lamb shift, by both electronic vacuum polarization and also 
by a hypothetical millicharged particle, can thus be approximated by taking the negative of the vacuum 
polarization energy shift of the (energetically lower) 25 level. When the resultant shift is divided by the 
observed discrepancy SE given in Eq. (jl.llbjl . one obtains the ratio 

f 00 dt 

X»= -V»(t)5p(t). (4.15) 



where, again with r = tj (4mV), 



a 1 



3tt SE 



n 



^) = ^(M^2S) = ^[l + ^|ll + ^l] . (4.16) 



a 3 mji 1 / 8v 2 t\ / 2w % /r x 



12tt SE V a 2 J V a 



Here, v = m e /mpi is ratio of the electron mass to the reduced mass of the muonic hydrogen system. The 
calculation of the ratio 

dt 



.■Va(t) 5 p{t) 

G(m M ) = — = ^joA (4.17) 

x " I jV,(t)Sp(t) 

then answers the following question: Suppose that the energy discrepancy SE in muonic hydrogen were due 
to the millicharged particle, then how much would the muon anomalous magnetic moment be changed by 
that same millicharged particle, in terms of the observed discrepancy 5a? Within the approximation (|4.8j) . 
the ratio G depends only on the mass (not on the charge) of the millicharged particle, 



oo 



dt , . 

jVa(t) 

m M )= ~;Z M dt , (4.18) 

because e 2 as given in Eq. (|4.8j) cancels. If we could find tum so that G{w,m) = li then a millicharged 
particle would be a serious candidate to explain both the muonic anomalous magnetic moment as well 
as the muonic Lamb shift discrepancy. 

The somewhat disappointing result of a numerical study of G(tom) is given in Fig. [3] In the mass range 
10 m e < mnj < 100 m e , the function G(mj\/) is in the range 25 < G(itim) < 300. The value G(itim) = 25 
implies that, if for given mass of the millicharged particle, the muonic hydrogen Lamb shift discrepancy 
is resolved, then we induce a discrepancy of theory and experiment for the muon anomalous moment by 
roughly 25 x 3.4 = 85 standard deviations. Expressed differently, the correction to the muon anomalous 
magnetic moment, expressed in units of the observed discrepancy <5a M , is larger by at least a factor 25 
than the modification of the muonic Lamb shift induced by the millicharged particle, expressed in units 
of the observed discrepancy SE. That means that a millicharged particle in the given mass range cannot 
explain both observed discrepancies. The observed discrepancy SE is too large to be explained by a 
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millicharged particle with m e <C tum "C "V because this would induce a prohibitively large modification 
of <5a M . Conversely, if a virtual millicharged particle provides an explanation for the observed discrepancy 
<5a M , then it will explain at most 4% of SE. Thus, a millicharged particle in the given mass range might 
still explain the discrepancy 8a lll but if that assumption is true, then the bulk of the explanation for SE 
has to come from a different effect (e.g., proton structure). 

We have previously stressed that muonic QED experiments are especially sensitive to a mass range 
m e -C rriM *C ni^ of the hypothetical particle. In principle, one might still explore the possibility of a 
millicharged particle with mass am e <C tum *C m e , because in that mass range, the correction to the 
hydrogen Lamb shift is still expressible in terms of a Dirac S function and therefore can be absorbed into 
a modified proton radius [7] . This study is indicated even if the mass range rriM *C m e is not the primary 
range tested by muonic experiments. According to Fig. 21 we find that, in principle, a particle with 
e 0.0179 and tum ~ 0.221 m e could explain both observed discrepancies 8E and Sa, simultaneously. 

However, this hypothetical particle is excluded for two reasons. First, a calculation of its contribution to 
the electron anomaly would be Sa e = 2.3 x 10 -10 which is much larger than the experimental uncertainty 
of ±2.8 x 10 -13 of the recent measurement [62] . A shift in the electron anomaly by Sa e = 2.3 x 10~ 10 would 
lead to a relative shift of the fine-structure constant by 2.0 x 10 -7 and thus, to a severe disagreement 
with other determinations of a (see Ref. [7]). The second reason is as follows. Because tum ~ 0.221 m e 
is below the electron mass, the effect of the virtual particle on the electronic hydrogen spectrum is no 
longer parametrically suppressed. A calculation shows that because the hypothetical virtual particle now 
is "too light," the proton radius inferred from the hydrogen spectroscopy experiments would increase to 
0.939(7) fm, because of the concomitant Dirac 5-likc vacuum polarization potential induced by the light 
millicharged particle. This is a priori not a problem, because the same modification (due to the mil- 
licharged particle) would have to be applied to the proton radius inferred from scattering experiments |S]. 
However, in that case, since the discrepancy SE is to be explained by the millicharged particle, the pro- 
ton radius inferred from the muonic hydrogen Lamb shift would be equal to the CODATA median value 
of r p = 0.8768 fm, and there would thus be a 9a deviation from the modified value r p = 0.939(7) fm 
inferred from electronic hydrogen spectroscopy. Therefore, the parameter range am e <C rriM "C rn e for 
the millicharged particle also can be excluded. 

4.3 Unstable Neutral Vector Bosons 

As evident from Figs. [3] and Fig. 21 our attempts to explain both muonic QED discrepancies with the 
simple form (|4.8[) of the modified vacuum polarization charge density have not proven successful. In prin- 
ciple, one may justify a more complicated ansatz for the modification, and with enough free parameters, 
it will certainly be possible to find a convenient representation that "explains" both discrepancies and 
is not in conflict with the electronic hydrogen Lamb shift and with the electron anomalous magnetic 
moment. This is not our goal. 

However, one further, specific form of a modified spectral density of the vacuum polarization deserves 
a discussion. A light, neutral vector boson has been investigated as a possible candidate to explain a 
prevailing discrepancy of the decay rate of orthopositronium (experiment versus theory, see Refs. |63L 
165]). which has eventually been resolved [66]. Just like the p(770) vector boson, such a hypothetical, 
virtual neutral vector boson would induce a small hump in the spectral density of vacuum polarization, 
corresponding to a resonance in the photon propagator. The modification would be restricted to a finite 
subintcrval of the t parameter. This possibility is not absolutely excluded by other searches [631164] because 
the vector boson might be unstable. We have performed extensive numerical experiments, for masses (and 
widths) of the virtual vector boson in the primary range m e <C my <C m p . One example is 

6p(t) = —Q(t- 20 m2) e ( 24m 2 _ t ) _ ( 4 19 ) 

For this choice, the muon anomalous magnetic moment discrepancy is reduced to 2.5 a, with the shifted 
theoretical prediction now lying above the experimental value. In addition, the muonic hydrogen Lamb 
shift discrepancy is reduced to 2.6tr, with the theoretical value still lying below the experimental one. The 
muonic hydrogen value of the proton radius would thus shift to a value of r p = 0.858 fm. The electron 
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anomaly a e is shifted by 4.2 x 10~ 14 which is below the current experimental uncertainty [55] ■ However, 
an evaluation of the additional vacuum polarization correction due to dp as given in Eq. (|4.19|) for elec- 
tronic hydrogen shows that the proton radius inferred from hydrogen spectroscopy would also have to be 
modified, namely to r p = 0.904(7) fm, because of the additional Dirac 5 potential induced by the virtual 
vector boson. This would leave a dissatisfactory 6.7a deviation between the modified radii from the two 
bound systems. Despite extensive numerical experiments, we have not found a simple, satisfactory pa- 
rameter combination that might explain both muonic QED discrepancies without significantly distorting 
the proton radius inferred from electronic hydrogen. 



5 Formation— Process Dependent Screening Corrections 

All hypothetical explanations for the discrepancy of theory and experiment in muonic hydrogen discussed 
so far in this article do not seem to lead to a satisfactory explanation. In order to understand a physical 
problem, null results can also be important, but are not gratifying and leave the effect unexplained. So, 
in order to fully understand the problem, we also analyze the experimental procedures used in Ref. [5]- 

In the experiment [5], devices have been installed in the beam line to extract electrons. E.g., as revealed 
in Fig. 2 of Ref. [5], the beam line is designed so that muons pass two stacks of thin carbon foils, and 
the electrons released from the foils are then extracted from the beam line via E x B drift and detected 
in scintillators. This leads to a separation of muons and electrons, and only muons continue in the beam 
line, to hit the gas target. The molecular hydrogen gas target is installed in the beam line behind the 
foils, i.e., after the electrons have been extracted from the beam (but not from the hydrogen molecules 
in the gas target). 

In many other experiments involving the high-precision spectroscopy of muonic transitions, the electron 
screening correction has been the limiting factor in analyzing the experiments [24(131] . Quite exotic 
processes involving muonic atoms have been studied in the literature (see, e.g., Ref. [67]). One of the most 
striking surprises is the role of exotic bound states like the well-known [68] molecular state composed of 
a pppT "nucleus" , and another proton, bound together by two orbiting electrons. The mentioned state is 
known to play a role in muon-catalyzcd fusion. One might thus ask to which extent the electrons in the 
H2 gas target may influence the observed lines. Note that the formation process of muonic hydrogen is 
complex and it is nontrivial to exclude the contribution of resonances from other muonic bound states; 
for this reason, resonances of the muonic molecules ppu. have been studied |69| . The first excited state of 
the pp/j, molecule is predicted to have a lifetime of 0.0713 ps, close to half that of a muonic hydrogen atom 
in the 2P state. As the width of the resonance observed in the experiment [5] is close to the calculated 
width for the muonic hydrogen 2P state, the authors of [5] exclude the possibility that the molecular 
resonance may contribute to the observed signal. 

However, these considerations do not exclude contributions from p[i~ e~ atoms composed of a proton, a 
negative muon and an electron, which might be formed in the gas target (these would be heavy analogues 
of the hydrid ion H~). Because the muon's orbit is close to the proton, the proton charge in pu~e~ is 
shielded from the outer electron. However, the muon and proton in the 2S or 2P state form a neutral 
core that interacts with the electron via dipole interactions. It is known that an ion-atom interaction with 
a functional dependence of the form — 1/r 4 can form bound states (for a recent numerical investigation, 
see Ref. [70]). Here, the "ion" is the electron, whereas the "atom" is the muon- proton core. The form of 
the interaction potential can be derived as follows. 

We first slightly generalize the problem and assume that the proton is a nucleus with charge number Z . If 
we have a nucleus with charge number Z and a muon and an electron bound to it, then the unperturbed 
Hamiltonian for the muon reads 

*„ = ^---, (5.D 
2m M r M 

and the unperturbed Hamiltonian for the electron is 

fle = jL_(^i)£. (52) 

2m R r R 
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because the outer electron merely sees the screened charge of the inner core, which is the nucleus of 
charge number Z plus the negative muon of negative charge, effectively reducing Z by one unit. For a 
molecular hydrogen gas target, we can set Z = 1 (approximately), and H e is approximately equal to the 
free Hamiltonian but will be corrected by a — 1/r 4 interaction, as detailed below. The total Hamiltonian, 
including the muon-electron repulsion, then is 



The perturbation is 



2m fJ , r„, 2m e r K 



a 



H I = H-H e -H fl = + — . (5.4) 



r, 



' in 



We expand this expression up to the dipole term and assume that r e > r M , and obtain the interaction 
Hamiltonian, 



a a r 



Hi « 1 h a^r Xe-Xu, (5.5) 

r e r e rj 

where x e = r e /r e , and x^ — f^/r^. The first two terms cancel, and the third gives the dipole interaction. 
We start from unperturbed states with the electron in a state with quantum numbers \n e £ e m e ) and 
the muon in the ^S)^ state. A calculation in second-order perturbation theory then gives the energy 
perturbation 



SE = - — (n e £ 




\{2S\rl\n»Pm ll ) \ 



2 

li 



P - ^2S 



(5.6) 



Here, the superscripts and subscripts identify the particles (e stands for the electron, and fi stands for 
the muon). The |2S') A1 state undergoes virtual transitions to muonic \n^P) states, and the term in curly 
brackets in Eq. (|5.6j) is recognized as the static polarizability of the 2S state. The energetically closest 
state to the muonic 25 state is the \2P 1 / 2 ) t j, state, which is energetically removed from 125")^ only by the 
Lamb shift. The energy difference is of the order 

4i /2 -4^^m,- (5.7) 
The dipole matrix elements for the virtual transitions of the muon are of the order of 

(2S\ri\n ft Pm fl ) — . (5.8) 

We assume that the electron is in a (superposition of) states with quantum numbers \n e £ e m e ) whose 
dimensions are of the order of the (electronic) Bohr radius, and so 



n, 



The result of our order-of-magnitude estimate is 



n e 4 ) ~ K) 4 • (5.9) 



SE - -a 2 (am e ) 4 ( -g— ^ ) = -0.42 meV (5.10) 



for the energy shift of the 2S state, if an additional electron is present. As this is the lower state of the 
2P-2S Lamb shift transition, the negative of SE, i.e., +0.42meV, needs to be added to the transition 
frequency, potentially explaining the discrepancy. The corresponding effect, evaluated for the muonic 
2P-L/2 state, has the opposite sign and therefore shifts the transition in the same direction as the effect 
calculated here. 

The main result of the above order-of-magnitude estimate is as follows. If for some reason, an electron is 
spatially displaced from the muonic hydrogen atom by only a few Bohr radii at the time of the laser-induced 
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Lamb shift transition, then the static polarizability of the 2S muonic state induces a systematic shift of the 
muonic transition that has the right sign and might be large enough to explain the observed discrepancy. 
In order to assess the viability of the pfi~ e~ atom hypothesis, one would have to calculate its spectrum, 
its ionization cross sections in collisions in collisions with other hydrogen molecules in the gas target. 
Furthermore, it would be necessary to study its inner Auger rates of p[i~ e~ as a function of the state 
of the outer electron, and its production cross sections in the collisions that take place in the molecular 
hydrogen target used in the experiment This is beyond the scope of the current article. The occupation 
numbers may depend on the formation process. In order to explain the single, well defined resonance 
line seen in the experiment JSj/, one would have to assume that formation proceeds predominantly into a 
specific state of the p/jTe" atom (or into a sufficiently narrow subset of resonances) . These requirements 
may, in the end, reveal that the current hypoethesis cannot explain the observed shift of the resonance 
line from the predictions of QED theory. Thus, we do not claim that the calculated effect necessarily 
explains the discrepancy of theory and experiment. However, as we were unable to discern viable theoretical 
explanations, our general statement is that it may be worthwhile to study possible explanations based on 
a process- dependent effect. 

Another mechanism by which an electron screening correction could enter the analysis of the experi- 
ment [5] might be from neighboring hydrogen molecules or atoms. In general, one assumes [7T] that the 
initial capture into highly excited states with principal quantum number n m *J m^/m e m 14 takes place 
via the reaction fi~ + H2 — > (ppfj,~e~) + e~ , where the muon is captured into a molecular orbit and an 
electron is ejected. Auger rates induced by neighboring hydrogen atoms and molecules dominate in the 
deexcitation process for liquid hydrogen [72], but are suppressed for a less dense hydrogen gas target. 
The deexcitation and muonic hydrogen formation process has been analyzed both theoretically [71H75] 
as well as experimentally [761177] . Still, there is a multitude of available reaction channels and conceivable 
bound states. It may thus be worthwhile to reexamine the dynamics of the formation process of muonic 
hydrogen in the 5 T magnetic field with a special emphasis on the hypothetical existence of neighboring, 
screening electrons that may perturb the observed frequency in the experiment [5]. 

6 Conclusions 

Recently, two serious discrepancies of theory and experiments have been observed for muonic QED sys- 
tems: the muon g factor discrepancy of 3.4 a appears to persist [H3], and for the muonic hydrogen Lamb 
shift [5], an even larger discrepancy is observed which amounts to 5.0 a. The proton radius inferred from 
the muonic hydrogen Lamb shift is 5.0 a smaller than that determined by electronic hydrogen spec- 
troscopy [7] , or alternatively, if the proton radius from hydrogen spectroscopy or from electron scattering 
is used in order to obtain a theoretical prediction for muonic hydrogen, we observe a 5.0ct discrepancy 
of theory and experiment. The muonic Lamb shift discrepancy is statistically more significant and thus, 
arguably, more urgent to be resolved than the muon g factor. Proton structure effects (discussed in 
Sec. [3]) and small modifications of the vacuum polarization potential (see Sec. 2]) might be discussed as 
hypothetical explanations for the observed discrepancy. 

As shown in Sec. El an elaborate reevaluation of the two-photon exchange graph for muonic hydrogen, 
including its inelastic (proton polarizability) contribution, might contribute to an explanation of the 
muonic Lamb shift discrepancy. The inelastic part has been evaluated [TT]|36l|37] *° ^ e on the order of 
0.0015(4) meV. Unless a reevaluation reveals a somewhat surprising enhancement of the contribution by 
at least an order of magnitude, the discrepancy of 6E = 0.31 meV will persist. An alternative explanation 
due to a conceivable "dip" in the slope of the proton form factor in the momentum range probed via 
muonic hydrogen spectroscopy (see Sec. I3.1j) could only be excluded conclusively by a measurement of 
the proton form factor in scattering experiments probing the momentum range indicated in Eq. (]3.3p . 
While certain dip and hump structures in the proton form factor have been seen in experiments [5] and 
in theoretical calculations [78], one has to admit that the "dip" hypothesis seems somewhat remote. 

In Sec. El we investigate the role of a hypothetical millicharged particle in a numerically small, but 
important modification of the spectral density of vacuum polarization. We focus on millicharged particles 
and do not treat hypothetical supersymmetric models. The conclusion is as follows: The most immediate 
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theoretical ansatz for the modification 5p(t) given in Eq. (|4.8[) comes from a low-energy millicharged 
particle; a model for a hypothetical virtual natural vector boson resonance is described by Eq. (|4.19j) . Our 
numerical experiments show that neither a simple modification of the vacuum polarization function due to 
a millicharged particle nor due to an unstable intermediate vector boson can simultaneously explain both 
discrepancies observed for the muon anomalous magnetic moment and for the muonic hydrogen Lamb 
shift without significantly distorting the proton radius inferred from electronic hydrogen. Millicharged 
particles and virtual vector bosons could explain the muon anomalous magnetic discrepancy while having 
a negligible effect on the muonic Lamb shift, but not vice versa. If these particles are unstable, then other 
parametric bounds (e.g., those following from the experimental investigation reported in Ref. |48j ) may 
not be applicable to an analysis of their virtual contributions within vacuum polarization loops. 

Our investigations severely restrict the parameter space available for a modification of the vacuum po- 
larization spectral function due to either millicharged particles or unstable virtual vector bosons in the 
low-energy domain. Being statistically (barely) significant, the muon anomalous magnetic moment dis- 
crepancy is numerically so small that it restricts possible modifications of the vacuum polarization spectral 
function due to hypothetical millicharged particles to such small coupling strengths that a simultaneous 
explanation of the comparatively large observed discrepancy in muonic hydrogen becomes impossible. In 
a small parameter range, a virtual vector boson resonance modification somewhat reduces the statistical 
significance of the combined discrepancies in muonic hydrogen and for the muon anomalous magnetic 
moment, but at the cost of inducing a prohibitively large modification of the (electronic) hydrogen Lamb 
shift. 

Hitherto undetected virtual particles from extensions of the standard model might still explain the two 
discrepancies SE and Sa^, but a more complex structure of the modification of the vacuum polariza- 
tion spectral function would have to be generated, with — possibly — both attractive as well as repulsive 
modifications. Repulsive modifications of the vacuum polarization potential have been discussed in the 
literature I601179|. A more complex modification of the vacuum polarization has more free parameters. 
Therefore, a conceivable determination of these parameters becomes possible, if at all, only when more 
experimental data on muonic systems (e.g., the muonic helium ion) become available. Our discussion in 
Sees. 14.21 and 14.31 does not exhaust all possible theoretical explanations from extensions of the standard 
model but covers particular models that "suggest themselves." A conceivable explanation from an exten- 
sion of the standard model is constrained by the comparatively large discrepancy in the muonic hydrogen 
Lamb shift, the comparatively small discrepancy in the muon g factor, and must respect the excellent 
agreement of theory and experiment for the electron g factor and the hydrogen Lamb shift. 

Having discussed the most attractive and far-reaching theoretical consequences of the recently observed 
discrepancy [5] and having found negative results, we then proceed to indicate a candidate effect for 
a process-dependent correction to the energy levels that may be relevant to the experiment [5] (see 
Sec. [5]). We believe that a study of exotic bound states and a reexamination of the dynamics of the 
formation process indicated in Sec.[5]might be interesting in its own right even if eventually, more detailed 
calculations might show that these effects do not offer the explanation for the observed discrepancy in 
muonic hydrogen. 

High-precision QED experiments arc excellent probes of conceivable low-energy modifications of the 
standard model. Our considerations highlight the need for further experimental evidence regarding muonic 
bound systems, before conclusive statements are possible. As discussed in Sec.[H discrepancies of theory 
and experiment in muonic systems have been encountered in the past and have eventually been resolved. 

Even if theoretical considerations and careful considerations of additional systematic effects concerning 
the experiment [5] fail to resolve the proton radius discrepancy, then there is a third way to resolve it, 
based on spectroscopic methods. It works as follows. In general, the fundamental constants derived from 
atomic spectroscopy are highly intertwined (not only the nuclear radii). One example is the Rydbcrg 
constant. If the new proton radius from Ref. [5] is inserted into the evaluation of electronic hydrogen 
spectra, then a very precise Rydberg constant is obtained, which is given as [5] 

c-Roo =3 289 841 960 251(5) kHz [1.5 x 10~ 12 ] , (6.1) 

where we indicate the relative uncertainty in square brackets. In an unpublished PhD thesis |80l completed 
at MIT (Cambridge) in 2002, based on Rydbcrg transitions which are manifestly independent of the 
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proton charge radius, the value 

cRoo = 3 289 841 960 306(69) kHz [2.1 x 10~ n ] , (6.2) 

is obtained [the difference to the value in Eq. (|6.1[) is 0.8u]. The CODATA value is 

= 3 289 841 960 361(22) kHz [6.6 x 10~ 12 ] , (6.3) 

and the difference to (|6.2[) also is 0.8<r, but "to the other side." The Rydberg state value (|6.2[) lies in 
between the values given in Eqs. (|6.1[) and the CODATA value (|6.3|) . Quite surprisingly, an improved 
measurement of the Rydberg constant based on ionic Rydberg states, though in itself independent of 
the proton radius, could thus settle the proton radius question connected with the muonic hydrogen 
measurement. Such a measurement is currently being pursued at the National Institute of Standards and 
Technology [ST] . 
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